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ABSTRACT 



Simple models of gravitational arcs are crucial for simulating large samples of these objects with full control of the input parameters. 
These models also provide approximate and automated estimates of the shape and structure of the arcs, which are necessary for 
detecting and characterizing these objects on massive wide-area imaging surveys. We here present and explore the ArcEllipse, a 
simple prescription for creating objects with a shape similar to gravitational arcs. We also present PaintArcs, which is a code that 
couples this geometrical form with a brightness distribution and adds the resulting object to images. Finally, we introduce ArcFitting, 
which is a tool that fits ArcEllipses to images of real gravitational arcs. We validate this fitting technique using simulated arcs and apply 
it to CFHTLS and HST images of tangential arcs around clusters of galaxies. Our simple ArcEllipse model for the arc, associated to 
a Sersic profile for the source, recovers the total signal in real images typically within 10%-30%. The ArcEllipse+Sersic models also 
automatically recover visual estimates of length-to-width ratios of real arcs. Residual maps between data and model images reveal 
the incidence of arc substructure. They may thus be used as a diagnostic for arcs formed by the merging of multiple images. The 
incidence of these substructures is the main factor that prevents ArcEllipse models from accurately describing real lensed systems. 

Key words, gravitational lensing: strong - methods: analytical 



1. Introduction 



Gravitational lensing is an important tool for a variety of as- 
trophysical applications. In particular, gravitationally lensed 
arcs produced by massive galaxy clusters can p rovide infor- 
mation on their total inner mass distribution ([Hattori et al. 
1999; Guzik & Seli aid 12001 iMandelbaum & Seliakll2006t iTreu 
201Q| iKneib & Nataraianl 1201 ll) . Statistical analyses of grav- 
itational arcs provide constraints on the evolution of sttuc- 
ture with cosmic time and on the cosmological parameters 
(iHilbert et alJl2007t IVuissoz et al J 120071: iKochanek et aljr 2006 ; 
i Meneghetti et al.l 120041: iGolse et all 120021: iBartelmann et al. 
1998: 'Estrada et al.l 120071: iHennawi et all l2008t iKausch et al. 
2010; Mo re et aD2012l) . Because of the flux amplification asso- 
ciated to gravitational lensing, these arcs are also a probe for the 
high-redshift universe, which is unaccessible through unlensed 
sources. 

On the other hand, gravitational arcs are very rare objects, 
since each square degree of the sky contains approximately 
one galaxy cluste r massive enough to p roduce arcs by grav- 
itational lensing (Glad ders & Yed 120051) . So far an order of 
hundreds of these objects were detected (jLuppino etal. 1999; 
Zaritskv & Gonzale z! 120031: Gladders et al.l 120031: ISand et all 
2005HCabanac et alj2007l:lBolton et al.l2008l) . Future wide-field 



Survey (DESfl will provide large samples of gravitational arc 
systems that will enable numerous statistical studies. Because 
of the large area, it will be virtually impossible to visually in- 
spect the whole survey to look for gravitationally lensed sources. 
Therefore, automated detection methods are ultimately needed. 

There are a few proposals for automated detection al- 
gorithms in the literature, although non e of them i s con- 
sidered as the default algorithm (e.g., Lenzen etall 2004: 
i Horesh et ail 120051: ISeidel & Bartelmannl 120071: ICabanac et al] 
l2007fc iMore et all 120121: iBom et al.l 120121) . Moreover, most of 
these algorithms are not fully characterized in the sense that 
they do not have a determined efficiency and contamination. 
These quantities could be quantified by using samples of im- 
ages with simulated arcs, which would act as truth tables for this 
type of study. Furthermore, simulated arc images are useful to 
test in conu'olled studies arc measurement tools and image pre- 
processing algorithms used to enhance their detectability. 

Currently, the literature lacks sufficiently detailed works re- 
lated to the fast simulation and characterization of gravitational 
arcs. Arc simulations have mostly been based on ray-tracing 
techniques given a distribution of source and lens properties. The 
common practice is to measure the arc's length L, width W, and 
its distance to the center of the lens. F ew authors have t ried to 
constrain their shapes and orientation. Miralda-Escudel (Il993h 
introduced the idea of measuring L by fitting a circle arc through 



surveys with sub-arcsecond seeing, such as the Dark Energy 



http ://w w w. darkenergy survey, org/ 
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the arc center and its extremes. Bartelm ann & Wei"ssl dl994) have 
pioneered the idea of approximating gravitational arcs as an el- 
lipse whose major axis is equal to L (the so-called ellipse fit- 
ting). However, we are not aware of any preceding study that 
has introduced a surface brightness distribution to simulate and 
characterize gravitational arcs. 

Motivated by these questions, we developed a code called 
PaintArcs which draws geometric figures that mimic arcs and 
adds them to images. It currently uses the prescription of an 
ArcEllipse, which is an analytical expression for the deformation 
of an ellipse such that one of its main axis becomes a circle seg- 
ment. Only a handful of fully controllable parameters go into this 
simple prescription. We use the ArcEllipse in conjunction with 
a Sersic profile to take into account the surface brightness distri- 
bution of the simulated arc. We also present ArcFitting, which 
is a code that uses the ArcEllipse+Sersic model to estimate the 
morphological and structural parameters of real arcs. 

The prescriptions that we adopt here for the shape and 
structure of arcs are certainly a simplification of real gravita- 
tional arcs, which are often formed by the merging of individ- 
ual multiple images of a source. One such situation are the gi- 
ant arcs formed when a source is close to the cusp of a caus- 
tic. Particularly when real arcs are imaged in high resolution, 
the ArcEllipse shape with a Sersic profile will not accommodate 
their substructures. Our goal here is to exactly to quantify the 
observed discrepancies from a simple model and, as a result, as- 
sess which improvements are most promising and to what extent 
these simple models can be used in controlled experiments for 
automatic arc detection algorithms. 

This paper is organized as follows. In §|2] we present the 
ArcEllipse, including the addition of a surface brightness pro- 
file and asymmetry. In |3]we describe the PaintArcs code, ac- 
counting for image pixelation, and including seeing effects and 
Poisson noise in the simulated image. The ArcFitting methods 
are presented in §|4] and their application to fit real images are 
presented in §0 Finally, in §|6] we discuss the results and future 
work. 

2. ArcEllipse 

The ArcEllipse consists of a simple idea to analytically express 
arc shapes as the distortion of an ellipse such that one of its main 
axes is bent into an arc of a circle. The resulting geometrical 
figure can be used to create a surface brightness map that mimics 
a gravitational arc. 

An ordinary ellipse is formed by the set of points whose dis- 
tances to the x and y axes, weighted by the semi-axes, a and b, 
satisfy 



(1) 



Following the same idea, we can define the ArcEllipse as the 
set of points whose distances from a point on the circumference 
along the tangential direction (r c AO) and along the radial direc- 
tion (Ar) satisfy 



(r c A0 



a I \ b 



(2) 



where r c is the radius of curvature of the circle (and thus the 
curvature radius of the arc). 

This simple expression indeed gives a shape that resembles 
an arc, as we can see in the left panel of Fig. \T\ We define the 
length of an ArcEllipse as the length of the arc segment between 



the two extreme points where Ar = 0. This leads to L — 2a. The 
width at AO — is given by W = 2b. Therefore, the length-to- 
width ratio of the ArcEllipse is simply given as L/W = a/b. 

In polar coordinates, if the origin of the coordinate system 
corresponds to the center of curvature of the ArcEllipse, we have 
Ar = r - r c and AO — - 6q, where O is the position angle of the 
ArcEllipse center, such that 



r c (Oo-0) 



(3) 



The ArcEllipse is therefore defined by four parameters: the 
length 2a, the width 2b, the curvature radius r c , and the position 
angle 6 . 

From the solution above, it is simple to show that the area of 
the ArcEllipse is given by 



■*0o+a/r c rr + (8) 

AatcEII = | | rdrdO 



XVo+a/r c p 
j— fl/r c J r_ 

M,+«/r c (r + (6») 2 ) - (r_(6») 2 ) 

Je -a/r c 2 



dO ■■ 



= nab — 



nLW 



(4) 



i.e., identical to that of an ellipse with semi-axes a and b. 

In order to represent an arc, we must impose some limits on 
the ArcEllipse parameters such as b < r c , to ensure the arc shape, 
and a < 2nr c , to avoid closing the ArcEllipse into a ring. 

2.1. Surface brightness profile for ArcEllipses 

To create a surface brightness map with ArcEllipse isophotes, 
from a given radial profile, we need to simply replace the argu- 
ment in a radial profile by the function 



R = 



'r c (0-0 ) 



(5) 



which is constant over ArcEllipses. Notice that in this case the 
ArcEllipse must be centered at the origin (center of curvature). 
We can rewrite this expression as 



R=T ^[r c (0-0 )(l-e)] 2 + (r-r c ) 2 = yR 



prof » 



(6) 



where e = 1 — b/a is the ellipticity. 
2.2. Asymmetric ArcEllipse 

Since the gravitational arcs in general are not symmetric, we can 
also add some asymmetry to the ArcEllipse model. The simplest 
way to implement that is to consider two different values for the 
semi-major axis, a\ and a^, preserving the semi-minor axis as 
b; in other words, two ArcEllipses with different a matched at 
6q. Points whose angle are larger that 0q will therefore have 
ellipticity 



ei = 1 



whereas points whose angle are smaller that 6q will have 



e2 — 1 • 

a 2 



(7) 



(8) 
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°0 2 4 6 8 10 12 14 16 °0 2 4 6 8 10 12 14 16 

pixel pixel 

Fig. 1. Geometrical figures produced by the ArcEllipse. Left panel: symmetric ArcEllipse created with e = 1 - b/a = 
0.9, b — 1 pixel, r c - 15 pixels, and 6q - 45". Right panel: asymmetric ArcEllipse created with e\ = 0.9, e2 = 0.75, 
with the remaining parameters as in the left panel. 



Now the length is given by 
1 1 



L — a\ + fl2 — b 



1 -, 



1 - <?2 



(9) 



and Eq. © still holds. 

In the right panel of Fig. Q] we show an example of an asym- 
metric ArcEllipse. It has the same values of r c and Go as the sym- 
metric arc shown on the left panel. But now we set ei = 0.9 and 
e 2 = 0.75. 

3. PaintArcs: an implementation of ArcEllipse 

PaintArcs is a code written in Python that implements a model 
to simulate objects with arc morphology and adds them to as- 
tronomical images. It performs the digitalization of the surface 
brightness profile whose radial argument is given by the Eq. ©. 
It also performs point spread function (PSF) convolution, to add 
the effect of the atmosphere and instrument, and adds the effect 
of Poisson fluctuations on the image counts. The code makes 
use of functions from SLtoolfl. a library for image processing, 
catalog manipulation, and strong-lensing applications. 

PaintArcs is currently set to use the ArcEllipse model to 
mimic gravitational arc shapes. But it is built in a modular way 
that can accommodate other prescript ions. The sur face bright- 
ness distribution currently uses the ISersicI (119681) model, for 
which the intensity profile is given by 



I s (r) = I e exp\-b n 



(10) 



where I e is the intensity at the effective radius r e that encloses 
half of the total light from the model. The constant b„ is defined 
in terms of the parameter n, which describes the shape of the 
profile. The value of b„ can be computed numerically or obtained 
using analytical expressions that approximate its value. In this 
work we use the approximation of ICiotti & Bertinl (119991) . given 

by 



b„(n) 



1 4 
2n - - + 

3 405n 



46 



131 



255 15« 2 1148175« 3 



The SLtools library is available at http://che.cbpf.br/sltools/ 



2194697 



+ 0{n 5 ), 



306907 17750n 4 



which is better than one part in « 10~ 4 for n > 0.36. Note that 
if r e is defined in a different way, the expression for b„ will be 
different. 

For an arc created with ArcEllipse we choose b = r e . In 
other words, we distort the circular profile only in the tangential 
direction, preserving its effective size (r e ) in the radial direction. 
In this case, the ArcEllipse with Seisic profile is given by 



I(r,6) = Is(R pm f(r,6)) 



(ID 



= I Q exp\-b„ 



V[r e (0-0 o )(l-e)] 2 + (r-r c ) 2 



l/n\ 



where Iq = I e e h " . 

We need to normalize the intensity of the ArcEllipse I(r, 6) 
such that the total signal of the arc, given by 



■In 



I(r' ',ffy 'dr' 'dff ', 



(12) 



is equal to the total signal of the arc in the image, which is given 

by 



1 total 



jQ-0.4(m-m_ T ,) 



(13) 



where m and m zpt are the total magnitude of the arc to be sim- 
ulated and the zero point of the magnitude scale of the image. 
Therefore, the constant 7 can be found as 



/o = 



1Q-0.4(m-m ;pf ) 



VM9'-9 )(l-c)] 2 +(r'-)>) 2 



1/H 



.(14) 



r'dr'dff 



For asymmetric ArcEllipses, the denominator of Eq. [14] 
should be replaced by a sum of two integrals, associated to the 
two values of ellipticity (ei and see 82.2I >. each one on either 
side of 6*o. 

To illustrate the method presented above, we show an arc 
created following the process described above in the left panel 
of Fig. 
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Fig. 2. Example of an ArcEllipse created by PaintArcs using a Sersic profile. The left panel shows the pure arc, created using 
r c = 17", r e = 0.6", O = 40°, n = 1, e t = 0.85, e 2 = 0.75, m = 21, and m zpt = 31.83. The pixel scale of the image is 0.154". The 
middle panel shows the application of a Gaussian convolution (FWHM = 1.1") on the pure arc image. The right panel shows the 
effect of adding Poisson noise to the convolved arc. 



3.1. Pixelation 



3.2. Adding the seeing effect 



To create a digitized image from a continuous intensity model, 
we need to integrate the object signal over a discrete number of 
pixels. We thus need to control the change in signal caused by 
the pixelation and try to keep it lower than some tolerance value, 
e. In a typical position in our image, close to the Sersic effective 
radius r e , the profile intensity is expected to vary by 



or 



(15) 



where 6 is the linear size of the pixel. Therefore, the fractional 
error in the pixel intensity is 



dl 



dl 



I I(r=r c +r e ,S=6 ) & r 



(16) 



We keep this error within a tolerance limit by estimating the 
intensity at a number of subpixels. We can find the number of 
subpixels n p by imposing dl/I = e as follows 



1 



1 



di 



e/(r=r c +r e; e=fl ) & r 



r=r c +r e ,9=8 



(17) 



The number n p represents the minimum integer number of sub- 
divisions of a pixel required to make the error on the pixelation 
less than e. 

In the example of Fig. [2 we used e = 0. 1 as the upper limit 
for the error on the digitalization of the arc signal. In this case, 
n p = 4, which means that each pixel was subdivided into 16 
parts. The exact value of the digitalization error, computed after- 
wards by comparing the total signal to be simulated (S to tai, see 
Eq. [T3l with the signal that was effectively distributed, is on the 
order of 10~ 3 . This value is much lower than the adopted e, at- 
testing that the criterium given by Eq. (TTTt is based on a profile 
region where the variation in the arc intensity is stronger than 
average, while the exact value of the error is computed using all 
pixels. Therefore the actual pixelation error is much smaller than 
e. 



To add the seeing effect to a simulated image created with 
PaintArcs, we convolve it with a 2D normalized Gaussian PSF, 
given by 



PSF = g(x,y) = 



expj- 




} 


J/ ex p{- 


x 2 , y 2 
2a 2 , 2tr 2 


> dxdy 



(18) 



where <x v and cr y are the standard deviation in x and y directions. 
The option o~ x + o~ y allows an elliptical Gaussian. 

Assuming a circular Gaussian (cr = o~ x = cr v ), the standard 
deviation can be obtained from the FWHM of an actual image 
through 



FWHM = 2 V21n2o-. 



(19) 



We implemented a Python module to convolve images with 
a Gaussian kernel whose size is a multiple of the standard devi- 
ation («o-). 

In the middle panel of Fig. [2] we show the result of the 
Gaussian convolution on the arc image previously presented (arc 
in the left panel of this figure), where we use a seeing FWHM of 
1.1", which corresponds to cr — 0.47", and n a = 4. 

Attesting the code's flexibility and modularity, new imple- 
mentations of PaintArcs are being developed to use alternative 
seeing kernels, such as a 2D Moffat function. 

3.3. Adding Poisson noise 

In a quantum detector such as a CCD the photon counts follow 
a Poisson distribution, 



P{x,ti)=n x —, 



(20) 



where fi is the mean of the distribution. 

We add Poisson noise to an ArcEllipse produced by 
PaintArcs by randomly picking values from the distribution 
given by equation d20b with the mean equal to the signal in each 
pixel of the simulated image after convolution with the PSF. This 
is done with a separate module that generates random numbers 
from a Poisson distribution. 
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In the right panel of Fig. [2] we show the result of adding 
Poisson noise to the arc of the middle panel of the same figure. 

3.4. Applying PaintArcs 

After the pixelation of the surface brightness profile of the arc 
created with PaintArcs with a Sersic profile and adding seeing 
and Poisson noise to its image, PaintArcs can paste this simu- 
lated arc onto a real image (for example one containing a galaxy 
cluster), which can be specified in its input configuration file. In 
this case, PaintArcs adds the arc model parameters to the image 
header for bookkeeping purposes. 

Figure [3] shows an example of an arc simulated with 
PaintArcs. It was added to the image of the galaxy 
cluster SOGRAS0328+0044, which was i maged by the 
SOAR Gravitational Arc Survey (SOGRAS, iFurlanetto et all 
120 121) .Two real arc candidates are also found around the clus- 
ter's brightest galaxy (indicated by arrows), and may serve as a 
visual guideline for the accuracy of the PaintArcs outputs. The 
simulated ArcEllipse is circled in the figure and was generated 
with the parameters m = 21.5, e\ = 0.9, e 2 = 0.8, r c = 7.7", 6q = 
40°, n — 3, and r e = 0.8". It was convolved with a Gaussian PSF 
whose FWHM was obtained from the image (FWHM = 0.7"). 
The origin of the coordinate system of the arc coincides with 
the center of the cluster's brightest galaxy. The pixel size of this 
image is 0.154". 




Fig. 3. Example of the application of PaintArcs on a real 
image: the circled arc was created by the code using an 
ArcEllipse+Sersic model with m = 21 .5, e\ = 0.9, e 2 = 0.8, r c 
7.7". 8 = 40°, n = 3, r e = 0.8", FWHM = 0.7" and was added 
to the image of the galaxy cluster SOGRAS0328+0044, where 
two real arc candidates (indicated by arrows) were found. 



In Section [2] we showed that the ArcEllipse depends on the 
following parameters (see Eq. |6]l 

- xq: x position of the point where the arc is centered (pixel) 

- yo'. y position of the point where the arc is centered (pixel) 

- r c : distance of the arc center relative to (xo,yo) (pixel) 

- 6q: orientation of the arc center relative to the x axis (radians) 

- e\ and e 2 . ellipticities of the arc 

- b: minor semi-axes of the arc (pixel). 

If the surface brightness distribution of the arc follows a Sersic 
profile, we add three more parameters to the fitting: 

- r e : radius that encloses half of the total luminosity of the 
Sersic profile (pixel) 

- n: index that describes the shape of the Sersic profile 

- Iq: intensity of the Sersic profile at the center of the arc (r = 
0). 

If we choose b = r e , as we do in PaintArcs, we ob- 
tain nine fitting parameters. Therefore, the set of parameters 
that represents the ArcEllipse with a Sersic profile is p = 

(x Q ,y , r c , 6 , e u e 2 , r e , n, 7 ). 

4.1. Creating an arc image to be fitted 

The ArcFitting code works on a postage stamp image of an ob- 
served arc (i.e., an image where the pixels not belonging to the 
arc are set to zero), which we call the arc image. To create the 
arc image we first run SExtractor 2.8.6 (Be rlin & Arnoutsl l996) 
on the original image that contains the arc. SExtractor then cre- 
ates a segmentation image, which assigns a single value (cor- 
responding to the sequential number of the object in the output 
image catalog) to all pixels that it associates to the object. The 
arc image is then created by copying all arc pixels defined by the 
segmentation image, but using their values as in the original im- 
age. All other pixels in the arc image are assigned a zero value. 
The size of the arc image in the x and y directions are multiples 
of the arc's major and minor axis values as found by SExtractor. 
The procedure described above to create the arc image makes 
use of functions from SLtools. 

4.2. Parameter determination from the arc image 

In the following we describe the measurements performed by 
ArcFitting to find the ArcEllipse parameters. But the ArcFitting 
is not only restricted to the model parameters, which currently 
are the nine ArcEllipse+Sersic parameters listed above. It also 
attempts to measure other commonly used arc parameters, such 
as its length L and width W. 



4. ArcFitting 

Motivated by the need for tools for measuring arc parameters, 
such as its length, width, curvature, and surface brightness, we 
also developed a tool that fits these properties to real (or realisti- 
cally simulated) gravitational arcs. We named it ArcFitting. 

Currently, ArcFitting is based on the ArcEllipse model with 
a Sersic profile. Therefore, it performs measurements on actual 
arc images that lead to estimates of the ArcEllipse parameters. 
Similarly to PaintArcs, ArcFitting is also a Python-based code 
that was designed in a modular way to allow different shape- 
and profile models to be adopted in the future. 



4.2.1. Determining the ArcEllipse parameters 

First, ArcFitting finds the two farthest points of the arc image, 
p\ and pz. For that, it uses any reference point to find the first 
extreme as the farthest point from the reference one. It then finds 
the farthest point from this first extreme through the farthest-of- 
farthest method. The code then names the extreme point whose 
angle is larger (smaller) than 9o (see Eq.|22]below) as p\ (pi). 

After that, ArcFitting obtains the points along the bisectrix 
between p\ and pi- This line is perpendicular to the line segment 
that connects these two extreme points. It then finds p^, which is 
the mean point of the arc along the b direction. 

This process is repeated twice to obtain the points p\ and ps. 
P4 is the mean arc position along the bisectrix line between p\ 
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and pj, whereas p$ is the mean arc position along the bisectrix 
line between p2 a nd py This meth od is an application of the 
Mediatrix method dBom et al.ll2.012b for n\, = 2, where rib is the 
number of bisections. The length (L) of the arc is then computed 
as the sum of the line segments connecting the five points 



L = P1P4 + PAPi + PiPs + P5P2- 



(21) 



In the left panel of Fig. [4] we illustrate this process of trac- 
ing the arc image and determining its length. The five points 
Pi, i — 1, 5 are shown along the arc image. The image used as an 
example in the figure is taken from an HST arc sample that will 
be presented in §0 

The ArcEllipse parameters xq, yo and r c are found by deter- 
mining the unique circle that passes through the extreme points, 
pi and p2, and through the point with maximum intensity, p max . 
The latter point is found as the arc point with maximum intensity 
among those that lie within 0.5cr fe from the line segments of Eq. 
|2T1 where 07, is the dispersion of the points along the bisectrix 
line containing p^,. This procedure is intended to eliminate arc 
image pixels located far away from the traced arc, which may be 
associated to a cosmic ray, or to a residual from the process of 
object segmentation and deblending, as these pixels would bias 
the Pmax estimate. The determination of p max is depicted in the 
right panel of Fig. [4] The pixel of maximum intensity in the en- 
tire image, m.\, failed to satisfy the 0.5<X£ distance criterion, in 
this case measured from the P1P4 segment. The second maxi- 
mum, ni2, however, is well-placed along the observed arc. 

The width of the arc is taken as W — 4ot,. The ArcEllipse 
semi-minor axis b is then taken as half of the width (b = W/2 = 
2cr b ). 

The angle Oq is obtained using the expression 



arctan 



(y c - yo) 

(X c - Xq) 



(22) 



where (x c ,y c ) are the coordinates of the p max point. 

The asymmetric ArcEllipse ellipticities e\ and e2 are ob- 
tained using Eqs. [7] and [8] respectively, where a\ is the size of 
the line segment p \ p max and «2 is the size of the line segment 

PmaxPl • 

The procedure outlined above is different compared to previ- 
ous works. Some studies, particularly of observed arcs, quote 
L, L/W, or other arc parameters, without clearly mentioning 
how they are determined. Other sources use no more than three 
points to infer L, either using a sum of arc se gments , as in 
eq. |2Tj or fitting a circular arc through them (lOguri 20021 
Dala Tet alJl2004t lLi et al] 120051: ISand et al.ll2005b Horesh et al " 



2005; Meneghetti et al. 2008; Horesh et al .l 2010t More et al 
2012). iLenzen et all (120041) and lKausch et a l. (2010) measured 



momenta over arc images to derive their sizes. 

Width measurements are usually based on the area of 
the arc divided by its length assumi n g so me arc geometry, 
usually a rectangle. Meneghetti et al.l d2008l) . however, mea- 
sured widths at several different positions along the arc length 
and assumed a median value. Some authors also account 
for seeing effects when m easuring W (Menegh etti et al.l f2008; 
IZaritskv & Gonzalezll2003l) . 

Arc orientation a nd curvature are not as common to find. 
One such reference is Luppi no et al.l ([1999), who used the chord 
length from the arc extremes and its sagitta l depth to est imate 
r c . Estimates of 9q ca n also be found in iQguril (120021) and 
iMiralda-Escudel (Il993l) . 

We are unaware of previous studies that accommodate a pos- 
sible asymmetry in the arc's shape. Finally, another element 



lacking in previous attempts to measure arc properties is the sur- 
face brightness distribution. This is the topic of the next subsec- 
tion. 

4.2.2. Determining the Sersic parameters 

We use x 1 minimization to determine the Sersic profile parame- 
ters, Iq, n and r e . This^ 2 minimization works as follows. 

We consider a data set that represents the information of N 
pixels in an arc image. We denote Ii mg (j) as the count of the 
j ,h pixel of the arc image and crj as its associated uncertainty, 
where j = {l,N). We also denote I arc (j,p) as the theoretical 
count of the j th pixel from ArcEllipse with a Sersic profile. This 
theoretical intensity depends on the pixel position and on the 
set of parameters p. Except for /o, n and r e , whose values are 
unknown, the other parameters are fixed by the procedure de- 
scribed in Sect. l4.2.TI Assuming that the pixel counts h mg {j) are 
independent from each other, we must minimize the quantity 



X 



-Z 



{larcij, P) ~ limgijif 



0~ 



(23) 



to estimate the remaining parameters that provide the best fit to 
the data. 

In general, we use the reduced x 2 



2 X 

Xred = 



N-l 



(24) 



as a measure of the concordance between the model and the data. 
If this value is much higher than unity, the model is not consid- 
ered a good description of the data. 

We assume that the pixel counts follow a Poisson distribu- 
tion, such that 



and Eq. (l24b becomes 



2 

Xred 



1 ^(larc(j,P)-Iimg(j)f 



(25) 



(26) 



The xf e d function minimization is carried out with 
Pyminuifl which is an extension for Python modules of the 
MinuiQ functions. 

The initial guesses for Iq and r e are determined through mea- 
surements on the arc image. We consider the intensity of the 
point p„ mx as the initial estimate for /o- The semi-minor axis b, 
whose estimate was described above, is the initial guess for r e . 
Given the lack of a direct method to estimate the parameter n 
from the image, the initial guess of this parameters is found by 
using the algorith Scan from Pyminuit, which calculates x^- ed for 
a range of n values keeping all other parameters in the p set fixed, 
i.e., the morphological parameters as determined in the previous 
subsection, and Iq and r e fixed at their initial guesses. 

After determining the initial guesses for /o, n and r e , we use 
the Pyminuit algorithm Migrad, which is an optimized gradient- 
based minimum search algorithm, to perform the jfi , minimiza- 
tion and to find the best-fitting values for these parameters. 



3 http ://code . google . com/p/py minui t/ 

4 http://seal.web.cern.ch/seal/snapshot/work- 



packages/mathlibs/minuit/ 
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Fig. 4. Illustration of the process of tracing a real arc in its image. The left panel shows the five p points used to compute the 
arc length. The right panel shows the high-intensity m points (dashed circles) used to find p max for the same arc. The center 
and radius of curvature are determined as the circle containing the points p\,p2, and p max = mi- The width is based on the 
line through the central point p^, whereas the asymmetric major-axes a\ and «2 are the lengths of the segments connecting 
'«2 = Pmax to the arc extremes p\ and p2- 



The approach outlined here to estimate the shape parame- 
ters directly from the image and to apply a x 2 minimization 
to the Sersic profile proved to be the most effective. We tried 
other alternatives, such as x 2 minimization (using Migrad) of all 
shape- and Sersic parameters, or an iterative method that com- 
bines Migrad for some parameters and brute force (Scan) for 
others. In many of the cases we tried, these methods did not con- 
verge, or resulted in a final much larger than unity. These 
failures are probably caused by a combination of strong degen- 
eracies among the ArcEllipse and Sersic model parameters and 
by the limited degrees of freedom when considering the number 
of available pixels in the arc image. 

4.3. Testing the ArcFitting algorithm 

To test and validate the ArcFitting code, we applied it to images 
of pure arcs (noiseless and with no seeing convolution) produced 
with PaintArcs, using the ArcEllipse+Sersic model. 

As an example, we considered an arc with true parameters 
as shown in the middle column of Table Q] The values recovered 
by ArcFitting are displayed in the right column, showing an ex- 
cellent agreement with the input values used by PaintArcs. This 
situation is typical of the validation tests. Indeed, ArcFitting was 
very successful in recovering the ArcEllipse and Sersic parame- 
ters in this idealized case. This is expected since both PaintArcs 
and ArcFitting use the same analytical expression for the arc 
shape and surface brightness profile, but demonstrates the valid- 
ity of the whole process of determining the morphological pa- 
rameters and using the initial guesses for scanning and fitting 
the Sersic parameters. 

5. ArcFitting application to real and simulated arcs 

We applied the ArcFitting code to three different sets of arcs. 
The first one consists of realistically simulated arcs. The other 
two sets are of real arcs, one from ground-based images and the 
other from HST images. 

We used the residual signal contrast AS/5, defined as 

XUimgti)- Iarc(j,P)) 

— = . (27) 

S 2 limgU) 

j 



Table 1. Input and fitted ArcEllipse and Sersic parameters for an 
arc created with PaintArcs using the ArcEllipse prescription and 
a Sersic surface brightness profile. 



Parameter 


True value 


Recovered value 


xo (pix) 


100.00 


98.50 


yo (pix) 


100.00 


96.88 


r c (pix) 


50.00 


52.85 


r e (pix) 


5.20 


5.96 


0o (rad) 


0.70 


0.72 


I (counts) 


100.00 


102.81 


n 


1.00 


1.06 


e\ 


0.90 


0.88 


e 2 


0.80 


0.77 



and the reduced x 2 , given by Eq. (1261 1. as metrics for the con- 
cordance between the ArcEllipse+Seisic model, resulting from 
ArcFitting, and the data. Both quantities are computed using 
only the pixels that belong to the arc as described in §14.11 

5.1. Applying ArcFitting on arcs simulated with AddArcs 

For this first application we used a set of arcs simulated with 
the AddArcs code (Brandt et al., in prep.), which was originally 
designed for Dark Energy Survey image simulations. The code 
uses as basic inputs a catalog of lenses, a catalog of sources, 
the cosmological parameters, and the properties of the images 
on which the arcs will be simulated (e.g., pixel scale, magnitude 
zero points, etc.). For the simulations used in this work, the in- 
put lens catalogs were give n by the Carmen Las Damas N-body 
cosmological simulation^] (McBri de et alj|201lh and the lenses 
are modeled as a projected Navarro-Frenk- White (tNavarro et alj 
1996. 119971 hereafter NFW) profile with elliptical mass distri- 
bution. The redshift and mass of each lens are taken from the 
catalo g, but the NFW c oncentration parameter is taken from the 
fits of iNeto et al.l (120071) . since the simulations do not have suffi- 
cient precision to determine it accurately. The sources are mod- 
eled as elliptical Sersic brightness distributions whose parame- 
ters are obtained from fits to galaxies in the Hubble Ultra-Deep 
Field (UDF) given by ICoe et al.l (|2006), which provides the in- 

5 http://lss.phy.vanderbilt.edu/lasdamas/ 
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put source catalog. This same catalog contains the photometric 
redshift information from multi-waveband data on the UDF. 

The sources are divided into redshift bins whose spatial dis- 
tribution is uniform in each bin. For each lens (i.e., dark mat- 
ter halo) i n the catalog, a number of sources is lensed using the 
Gravlens (Kee tonll200lh software. The number of sources to be 
lensed is obtained from the surface density distribution of the in- 
put source catalog. A first selection of sources is made using the 
local magnification computed with Gravlens and then the same 
code is used to produce a full brightness distribution by lensing 
finite Sersic sources from the UDF catalog. Finally, images are 
identified and are dimensionally estimated by SExtractor, and 
only elongated images are selected by the code. We selected a 
subset of nine arcs generated from these simulations to have a 
similar sample size as the two other samples. 

We ran ArcFitting on the arcs produced with AddArcs to in- 
fer their ArcEllipse+Sersic parameters and recreate them with 
PaintArcs. In Fig. [5] we show some examples of this applica- 
tion of the ArcFitting code. In the left panels we show the orig- 
inal simulated arcs (noiseless and without seeing convolution) 
from AddArcs. In the middle panel we show the ArcEllipse de- 
rived from ArcFitting and reproduced with PaintArcs. In the 
right panel we show the residual image resulting from the dif- 
ference between the two previous ones. All panels use the same 
grayscale. Some substructures are visible in these latter images, 
indicating the limitations of the ArcEllipse prescription. But the 
residuals are much less intense than the fitted arcs. 

We computed the residual signal contrast AS /S, as given 
in Eq. [27] between each AddArcs simulation and its fitted 
ArcEllipse+Sersic. In the left panel of Fig. [6] we show the re- 
sulting distribution of AS/5 values. The median AS/S value 
is -0.09, indicating that the ArcFitting process based on the 
ArcEllipse+Sersic model typically reproduces 91% of the signal 
of the input arc image. In the right panel of Fig. [6]we present the 
distribution of xL d resulting from the application of ArcFitting 
on arcs simulated with AddArcs, whose median is . = 5.27. 






\ 








\ 




\ 

\ 



Fig. 7. Selected ArcFitting results for HST arcs. Left panels: 
original arc; middle panels: arc created with PaintArcs using the 
parameters derived from ArcFitting applied to the original arcs; 
right panels: difference between original and reproduced arcs. 
From top to bottom: arc C8 found in Abell 68; arc HO of Abell 
963; arc M3 of Abell 2218, and arc B0/B1/B4 of Abell 383. 



5.2. Applying ArcFitting on arcs observed with the HST 

We here test the reconstruction of HST arcs using the ArcFitting 
code based on the ArcEllipse+Se rsic model. A samp le of HST 
arcs was selected from the work by Smi th et al.l d2005). We chose 
bright and isolated arcs with spectroscopic redshift measure- 
ments that confirm their nature as lensed sources. A total of 12 
arcs were selected. We use the same notation here as the orig- 
inal authors. The arcs are C4, C8, and C12 from the Abell 68 
galaxy cluster; HO, HI, H2, and H3 from Abell 963; B0/B1/B4 
from Abell 383; Ml and M3 from Abell 2218; and P0 and P2 
from Abell 2219. Arc properties such as magnitude, L/W, and 
distanc e to the BCG of these arcs are presented in San d et aD 
(12001 . 

In Figure [7] we present some results of the application of 
ArcFitting to our selected HST arcs. The panels of Figure Qare 
similar to those of Fig. [5] The left panels show the original HST 
arc images, the middle panels show the fitted ArcEllipse+Sersic 
model arcs, which was reproduced with PaintArcs, and the right 
panels show residual images. 

From the middle panels we can see that ArcFitting is effi- 
cient in recovering arc curvature and orientation. In some cases, 
as for the second and third arc in the figure, apart from the effi- 
ciency in recovering the arc shape, ArcFitting is also efficient in 
reproducing the arc brightness distributions, so that the residual 
images contain only little substructure and noise. However, for 



arcs with many substructures, such as the first and the last of that 
figure, ArcFitting cannot effectively reproduce the arc brightness 
profiles. 

We again used the residual signal contrast AS /S , given by 
Eq. [27] and the x 2 red , given by Eq. [26j to assess the agreement 
between model and data. The distribution of AS /S values for all 
12 arcs in our HST sample is shown in the left panel of Fig. [8] 
while the distribution of x 2 red is shown in the right panel of the 
same figure. The median of the AS/S distribution is -0.13, but 
values as high as -0.32 are found. The median xl ed is 18.39, 
but the distribution has a long tails toward much higher values. 
These results clearly attest that there is a stronger discrepancy 
between model and HST data than for the simulated AddArcs 
images. For all arcs, AS is negative, which means that ArcFitting 
is overestimating the arc signal. High values in both distributions 
correspond to the arcs that present many substructures (example: 
last arc of Fig. [7] which has AS/S = -0.32 and x 2 red = 76.26). 

In Fig. [9] we show a comparison of the L/W values measured 
bv lSandetal.1 (l2005h to those obtained with ArcFitting. Since 
the arcs HI, H 2, and H3 f r om A bell 963 are considered as a 
single object in ISand et aTl (|2005), these arcs are not included 
in this comparison. No strong systematics is seen, which shows 
that ArcFitting provides a robust and automatic measure of L/W 
ratio in real gravitational arcs. 
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Fig. 5. ArcFitting results for AddArcs arcs. Left panels: original arc; middle panels: arc reproduced with PaintArcs using the pa- 
rameters provided by the application of ArcFitting on original arcs; right panels: difference between original and reproduced arcs. 
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Fig. 6. Distribution of residual signal contrast AS/S (left panel) and ofx 2 red (right panel) resulting from applying ArcFitting on arcs 
simulated with AddArcs. 



5.3. Applying ArcFitting on arcs observed from the ground 

We also ran the ArcFitting code to reconstruct arcs observed 
from the ground. We selected bright and isolated arcs in cluster 
and galaxy group lenses from the CFHTLS S trong Lensing 
Legac y Survey (SL2S), presented in the work by Cabana c et al.l 
(120071) . Following the notation of the original authors, the 
seven selected arcs that compose our CFHTLS sample are 
SL2SJ021416-050315, SL2SJ021807-051536, SL2SJ085446- 
012137, SL2SJ141912+532612, SL2SJ142258+512440, 
SL2SJ143001 +554647, and SL2SJ143 140+553323. 

The results of applying ArcFitting on three arcs of our 
CFHTLS sample are shown in Fig. [10] As in the previous com- 
parisons, the left panels show the original arc images, the middle 
panels show the fitted ArcEllipse+Sersic model reproduced with 
PaintArcs, and the right panels show residual images resulting 



from the difference between the previous two. We can see that 
ArcFitting is efficient in recovering arc shape and brightness dis- 
tributions. 



In the left panel of Fig. QT| we present the residual signal 
contrast distribution for the selected CFHTLS arc sample. As in 
the HST arcs, we can see that ArcFitting overestimates the arc 
signal, since AS < for all arcs of this sample. The median of 
the AS/S distribution is -0.11 and AS/S = -0.20 in the worst 
case. In the rigth panel of the same figure, we show the distribu- 
tion ofx 2 rt>d , whose median is 10.76. In contrast to the HST arcs, 
no high x*. tail is seen; the only outlier in the ^f? , distribu- 
tion corresponds to the arc SL2SJ085446-012137, the same arc 
that presents the highest value in the AS/S distribution. The im- 
provement in ArcFitting of CFHTLS arcs with respect to HST 
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Fig. 8. Residual signal contrast (AS/S) distribution (left panel) andx 2 red distribution (right panel) resulting from applying ArcFitting 
on our selected sample of arcs from the HST. 
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Fig. 11. Residual signal contrast (AS/S) distribution (left panel) andx 2 red distribution (right panel) resulting from applying ArcFitting 
on our selected sample of arcs from CFHTLS. 
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Fig. 9. Compariso n between the L/W measurements presented in 
Sand et all d2005l) and those obtained with ArcFitting. The solid 
line represents the identity function. 



arcs is likely caused by the smearing of the arc substructures 
caused by seeing effects on the ground-based set. 

6. Conclusions and future perspectives 

We have presented and explored in detail the ArcEllipse, which 
is a simple analytical model that approximately reproduces the 
shape of real gravitational arcs. ArcEllipses are ordinary ellipses 
whose major axes are bent along an arc of a circle and which 
therefore mimic the shape of tangential arcs. They are fully char- 
acterized by a limited number of parameters that quantify their 
curvature, length and width, degree of deformation, and ori- 
entation. We also devised a simple prescription to incorporate 
asymmetries to their shapes. This description of arc morphology 
is more comprehensive than previously found in the literature. 
We also incorporated a Sersic light profile distribution to the 
ArcEllipses in an attempt to explore a rough structural model 
for lensed galaxies. 

We also implemented PaintArcs, which is a tool that cur- 
rently generates simulated ArcEllipse images, and ArcFitting, 
which is a code that derives model parameters from images of 
real arcs. Both are Python codes built in a modular way, thus 
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Fig. 10. ArcFitting results for CFHTLS arcs. Left panels: 
original arc; middle panels: ArcEllipse+Sersic image repro- 
duced with PaintArcs using the parameters resulting from 
ArcFitting; right panels: difference between the previous two 
panels. From top to bottom: arcs SL2SJ141912+532612, 
SL2SJ143140+553323, and SL2SJ021807-051536. 



allowing flexibility. PaintArcs, for instance, incorporates the ef- 
fects of pixelation, atmospheric seeing, and Poisson noise to the 
simulated ArcEllipse images, where the latter two can be turned 
on and off. ArcFitting is capable of deriving the model parame- 
ters using different techniques that vary from simple and direct 
estimates applied on an image to statistical inference methods. 
Even though both codes currently adopt the ArcEllipse shapes 
and Sersic profiles for the arcs, other more complex models can 
be added to their capabilities in the future. 

We applied the ArcFitting code to different types of gravita- 
tional arcs. Realistic arc simulations based on ray-tracing tech- 
niques applied to physical models of lens and source, such as 
those generated with Gravlens, have their general features and 
total signal reasonably well fit by an ArcEllipse+Sersic model. 
The resulting model from ArcFitting typically amounts to more 
than 90% of the original signal from the arc. Typical reduced^ 2 
values are on the order of 5. 

For real gravitational arcs around galaxy clusters, either im- 
aged from the ground or using the HST, the total signals are re- 
produced within 10-20% for ground-based images and within 
10-30% for HST images. Comparing the L/W values inferred 
from ArcFitting to those published for HST data shows that 
ArcFitting successfully quantifies their axial ratio. We attribute 
the poorer fitting results from HST images to the presence of 
substructure in the arcs. This opens up the possibility of using 
ArcEllipse models to enhance these substructures, helping in 
identifying giant arcs formed from merging of multiple images. 
Extending this reasoning, the ArcEllipse shape coupled with a 
single Sersic profile might in fact be used as a discriminator 
to filter out contaminanting elongated galaxies from a sample 
of arc candidates. Alternatively, gravitational arcs formed from 
multiple merged images of the same source may be fit by more 
than one ArcEllipse+Sersic component to improve their struc- 
tural modeling. 



Attempts to apply a full x minimization fail in general, 
given the relatively small number of degrees of freedom in the 
limited arc images available. Our results also indicate that the 
ArcEllipse prescription is too simplistic to yield significant im- 
provement in the models based on formal statistical inference. 
We thus reached a compromise set-up in ArcFitting by which 
statistical inference is applied only to the structural Sersic pa- 
rameters, while the shape ArcEllipse parameters are based on 
measurements over the arc images. The algorithm was succes- 
fully validated in this way by fitting a pure ArcEllipse+Sersic 
object with high precision. 

The usefulness of the ArcEllipse+Sersic prescription should 
not be underestimated. ArcEllipses of varying magnitudes, 
lenghts, L/W, and radii of curvature can be speedily generated 
by PaintArcs and used in a number of arc studies, taking advan- 
tage of the full control we have on the input arc parameters. The 
Dark Energy Survey (DES) simulations have already included a 
sizable sample of arcs generated by PaintArcs by the DES strong 
lensing working group. Simulated objects using PaintArcs allow 
an efficient testing of arcfinders, which are codes meant to au- 
tomatically detect gravitational arcs (Furlanetto et al., in prep.). 
In a 5,000 sq. deg. survey such as DES, usage of automated ar- 
cfinders is unavoidable, since their large areas preclude visual 
image inspection. 

The ArcEllipse also allows an approximate initial character- 
ization of observed gravitational arcs. Even such rough morpho- 
logical and structural parameters of observed gravitational arcs 
are not commonly found in the literature. Our approach is thus 
providing both a qualitative and quantitative improvement on 
the amount of information extracted from observed gravitational 
arcs. The efficiency of ArcFitting is also an important aspect. For 
large samples of arcs that contain perhaps thousands of them, as 
in the case of the DES, an efficient first pass on arc parameters 
is a necessary tool associated to the arcfinders themselves. 
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